Quantum dissipative Rashba spin ratchets 
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We predict the possibility to generate a finite stationary spin current by applying an unbiased ac 
driving to a quasi-one-dimensional asymmetric periodic structure with Rashba spin-orbit interaction 
and strong dissipation. We show that under a finite coupling strength between the orbital degrees 
of freedom the electron dynamics at low temperatures exhibits a pure spin ratchet behavior, i.e. a 
finite spin current and the absence of charge transport in spatially asymmetric structures. It is also 
found that the equilibrium spin currents are not destroyed by the presence of strong dissipation. 
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An opportunity to induce a net stationary particle 
current by unbiased external forces applied to a quan- 
tum dissipative one-dimensional (ID) periodic structure 
is provided when the system does not possess a center of 
inversion in real space Then the particle transport 
occurs due to the ratchet effect and the device works as 
a Brownian motor 2]. In the deep quantum regime the 
charge ratchet effect can only be achieved when at least 
the two lowest Bloch bands contribute to transport 

Recently a new research field of condensed matter 
physics, spintronics, has emerged. One of its central is- 
sues is how to generate pure spin currents (SC) in para- 
magnetic systems due to only spin-orbit interactions and 
without applied magnetic fields. Rashba spin-orbit in- 
teraction (RSOI) [i[ represents one of the possible tools 
to reach this goal since the spin-orbit coupling strength 
can be externally controlled by a gate voltage. One way 
to get pure SC is due to the intrinsic spin Hall effect 
[B|, Q expected in a high-mobility two-dimensional semi- 
conductor systems with RSOI [7|]. Such pure SC were 
experimentally detected through the reciprocal spin-Hall 
effect in Ref. An alternative is to induce pure SC 
through absorption of polarized light H. T he generation 
of pure SC by coherent spin rectifiers [10( has been dis- 
cussed only recently for a finite size setup with RSOI. 
However, the presence of dissipation has not been con- 
sidered up to now. 

In this letter we address the challenging task of how 
to implement devices which can work both as Brown- 
ian charge and spin motors. Here a natural and also 
principle question for spintronics arises: Is it possible 
to switch a device working as a charge ratchet to a pure 
spin ratchet mode where the charge current (CC) is com- 
pletely blocked? As mentioned above, when in a dissipa- 
tive system without RSOI transport is restricted to only 
one Bloch band, the charge ratchet mechanism does not 
exist 0. Whether the same effect takes place in a dis- 
sipative system with RSOI is an open and non-trivial 
question. In fact, the Rashba Hamiltonian is not in- 
variant under reflection of a transport direction. Thus 
the Rashba Hamiltonian itself already has a built-in spa- 




FIG. 1: (Color online) A schematic picture of the isolated 
asymmetric periodic quasi-lD structure described by the 
Hamiltonian In the center of the quasi- ID wire the peri- 
odic potential is weaker and gets stronger closer to the edges. 
Thus the electron group velocity is higher in the central region 
and tails off away from the center. 



tial asymmetry which due to the spin-orbit coupling can 
be further mixed with the periodic potential symme- 
try/asymmetry. The presence of dissipation additionally 
increases the complexity of the problem because the in- 
fluence of a dissipative environment on the orbital motion 
changes through RSOI the spin dynamics. 

In this work we focus on the moderate-to-strong dissi- 
pation case and address how to implement a device which 
under influence of unbiased external ac-driving yields a 
finite stationary spin current and at the same time blocks 
the directed stationary charge transport. To concretize 
our idea of a Brownian spin motor we consider a dis- 
sipative periodic system with RSOI and show that the 
spin-orbit interaction alone is not enough to produce SC: 
The system must additionally lack the spatial symmetry 
and its orbital degrees of freedom must be coupled. 

The full Hamiltonian of our problem is Hf u \\(t) — 
H + H cx t(t) + -Hbath, where H is the Hamiltonian of the 
isolated periodic system, H ext (t) describes an external 
driving, and -f/bath is responsible for dissipative processes. 

The isolated quasi-lD periodic system is formed in a 
two-dimensional electron gas (2DEG) with RSOI using a 
periodic potential along the x-axis and a harmonic con- 
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finement along the z-axis: 



H= 



2m 



ft 2 k°, 



to 



(pxk z - a z k x )+U 7 (x, z), (!) 



where Uy(x,z) — U(x)(l + jz 2 /L 2 ), k is related to the 
momentum operator as p = Kk, luq is the harmonic con- 
finement strength, fc so the spin-orbit coupling strength, 
U(x) the periodic potential with the period L, and 7^0 
the orbit-orbit coupling strength. This isolated structure 
is sketched in Fig. Q] as it could be realized by appropriate 
gate evaporation techniques applied to 2DEGs formed in 
III-V compounds. 

The periodic structure is subject to an external ho- 
mogeneous time-dependent electric field, E(t) = E(t)e x . 
It can be experimentally implemented using for example 
linearly polarized light. This yields H CK t — eE(t)x, where 
e is the elementary charge. We use the time dependence 
eE(t) = F cos(£l(t — to)), which is unbiased. 

The system is also coupled to a thermal bath. We as- 
sume the transverse confinement to be strong enough so 
that the probabilities of the direct bath-excited transi- 
tions between the transverse modes are negligibly small. 
Thus the environment couples to the electronic degrees 
of freedom only through x. Furthermore, in the spirit 
of the Caldeira and Leggett model [H}, we consider a 
harmonic bath with bilinear system-bath coupling. 

The dynamical quantities of interest are the ratchet 
charge and spin currents Jc,s(t) given as the statistical 
average of the longitudinal charge and spin current op- 
erators, Jc,s(t) = Tr[Jc,s/5(£)], where p(t) is the reduced 
statistical operator of the system, that is the full one with 
the bath degrees of freedom traced out. The CC operator 
is Jc(t) — —edx/dt and for the SC operator we use the 
definition suggested in Ref. Js(t) = d(a z x) /dt. 

It is convenient to calculate the traces using the basis 
which diagonalizes both x and <r z , because this requires 
to determine only the diagonal elements of the reduced 



density matrix. As shown in Ref. 13], for a periodic 
system with RSOI the energy spectrum can be derived 
from the corresponding truly ID problem without RSOI. 
This leads to so-called Bloch sub-bands. The 2DEG is as- 
sumed to be sufficiently dilute to neglect the Pauli exclu- 
sion principle in the temperature range of our problem. 
The upper limit of this temperature range is considered 
to be low enough so that only the lowest Bloch sub-bands 
are populated. The basis which diagonalizes x and a z 
becomes in this case discrete. The total number of the 
Bloch sub-bands is equal to the product of the number, 
-/Vb, of the lowest Bloch bands from the corresponding 
truly ID problem without RSOI, the number, 7V t , of the 
lowest transverse modes and the number of spin states. 
In this work we shall use the model with Nb = 1, Nt = 2. 
The total number of the Bloch sub-bands in our problem 
is thus equal to four. Using Nb = 1 we also assume that 
the external field is weak enough and does not excite 



electrons to higher Bloch bands. The representation in 
terms of the eigen-states of x for a model with discrete 
rc-values is called discrete variable representation (DVR) 
0, EH ■ Let us call cr-DVR the representation in which 
both the coordinate and spin operators are diagonal. De- 
noting the cr-DVR basis states as {\a)} and eigen- values 
of x and a z in a state \a) by x a and <r a , respectively, the 
CC and SC are rewritten as Jc(t) = — e J2 a x aP a {i) and 
J s(t) = Y, a a a x aP a {t), where P a (t) = (a\p(t)\a) is the 
population of the cr-DVR state \a) at time t. 

We are interested in the long time limit Jq° s of the 
currents Jc,s{t)i averaged over the driving period 27r/57. 

The advantage of working in the er-DVR basis is that 
real-time path integral techniques can be used to trace 

El El. More- 



out exactly the bath degrees of freedom 
over, at driving frequencies larger than the ones char- 
acterizing the internal dynamics of the quasi-lD system 
coupled to the bath, the averaged populations P a (t) can 
be found from the master equation, 
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T 0a P a {t), (2) 



valid at long times. In Eq. @ T a p is an averaged tran- 
sition rate from the state \/3) to the state a). 

The first task is thus to identify the cr-DVR basis. The 
eigen-states \l,kB,j, a) of cr 2 were found in [l]| for the 
case 7 = 0. The results obtained in [3] are straight- 
forwardly generalized to our model since for N t = 2 the 
operator z 2 (and any even power of i) is effectively di- 
agonal. The quantum numbers I, fee, j, o stand for the 
Bloch band index, quasi-momentum, transverse mode in- 
dex and z-projection of the spin, respectively. As men- 
tioned above I = 1, j = 0, 1. One further finds 



(l',k' B ,j',a'\x\l,k B ,j,a) = 

= Sj' t j6 a ' ta k' B + ak so \x\l, fee + cfc so )j, 



(3) 



where the index j under the bra- and ket-symbols in- 
dicates that the corresponding electronic states are ob- 
tained using the periodic potential U 7 j(x) = U(x)[X + 
"fh(j + l/2)/mu!oL 2 ]. For a fixed value of j the di- 
agonal blocks in Eq. © are unitary equivalent and 
thus the eigen- values of x do not depend on a. The 
eigen-values of the matrix k B \x\l, are analyti- 
cally found and have the form x^ >m j = mL + d^j, where 
to = 0, ±1, ±2 . . ., C = 1, 2, . . . , N B and the eigen-values 
d^.j are distributed within one elementary cell. Thus one 
can label the eigen-states of x as |£, m,j, a). The corre- 



sponding eigen-values are x^^ 



We see that 



the cr-DVR basis states \a) introduced above are just the 
|C, to, j, cr) states, that is {\a)} = {\^,m,j,a)}. 

To calculate CC and SC we use the tight-binding 
approximation assuming that the matrix elements 
(£', to', j', a'\H\Q, to, j, cr) with \m' — m\ > 1 are neg- 
ligibly small. Let us introduce the definitions for the 
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states \m,£) = |C = where {£} — {(j, a)} and 

£ = 1 O (0,1), £ = 2 (0,-1), € = 3 ^ (1,1), 
£ = 4 (1,-1). Correspondingly, we introduce hop- 
ping matrix elements = (m',£'|-ff]m, £) (m! m 

and/or ^ £) and on-site energies = (m,£\H\m,£). 

Due to the harmonic confinement and RSOI the sys- 
tem is split into two channels: one with £ = 1,4 and 
another with £ = 2, 3. The two channels are indepen- 
dent of each other, that is, transitions between them are 
forbidden. This picture is general and valid for an ar- 
bitrary number of the transverse modes. For clarity we 
below only consider the channel with £ = 1,4. Two in- 
dependent channels were also found for a different type 
of confinement in Ref. 

Assuming that the hopping matrix elements are small 
enough we can use the second-order approximation [3] 
for the averaged transition rates in Eq. ([2| . We have 



T$p= l / rfrexp 



h 2 



Q[t,J{u)]- 



m sm (-J 



(4) 



where x 



with 



dij. In 



Eq. (fj| J (x) denotes the zero-order Bessel function and 
Q[t, J (to)] is the twice integrated bath correlation func- 
tion being a function of time r and a functional of the 
bath spectral density J (to) 0,[l^. The dependence of the 
transition rates on the orbit-orbit coupling 7 comes from 
two sources. The first one is the Bloch amplitudes and 
the second is the difference Ad = d\$ — d\ : \. In a tight- 
binding model the periodic potential is strong and thus 
Ad can be made less than all the relevant length scales, 
Ad/l r <C 1, where Z r = min[L, y/h/muiQ, Ml/F]. Hence 

the main effect of the orbit-orbit coupling on T™ ^ m comes 
only through the Bloch amplitudes, and we neglect terms 
of order 0(Ad/l T ). 

We then arrive at the main results of our work, the 
absence of the charge transport, Jg° = 0, and the ex- 
pression for the non-equilibrium spin current (NESC), 
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r(0) 



AO) 



lulu 
sin[(fc B - k' B )L\lTo[F k , 



to 



in 



(5) 



where I? ^ , ij;? ^ are the integrals from (j4]) with and with- 
out driving, F ^ and F = 0, respectively, and 



x[<Y^ fcB - fe80 (rfi,i)<^ )fe ' +feTO (^,o)]*, 



(6) 



where ko (dij) is the DVR Bloch amplitude of the 

first band for electrons in the periodic potential U 7 j(x). 
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FIG. 2: (Color online) Non-equilibrium spin current, J^g, 
as a function of the amplitude, F, of the driving force for 
different values of the viscosity coefficient 77. Temperature 
fcBoitz.T = 0.5, spin-orbit coupling strength k so L — it/2, orbit- 
orbit coupling strength 7 = 0.1, driving frequency fi = 1. The 
inset displays the shape of the periodic potential. 



In Eq. ([5]) we have eliminated from Jg° the equilib- 
rium spin current (ESC), JJg, following Ref. [18]. The 
fact that the ESC turns out to be finite shows that the 
definition of SC suggested in Ref. [l2| does not automati- 
cally eliminate the presence of ESC. However, as pointed 
out in Ref. fl2| . this current really vanishes in insula- 
tors. This can be seen from Eq. When the potential 
is strong, electrons are localized, the dependence of the 
function Fk Bt k B on the quasi- momentum disappears, and 
as a result both ESC and NESC are equal to zero. This 
reasonable result is ensured by the spin current defini- 
tion taking proper care of the spin torque. It is inter- 
esting to note that ESCs are present even in a system 
with strong dissipation. As recently proposed in Ref. 
[l{|, ESCs can effectively be measured using a Rashba 
medium deposited on a flexible substrate playing a role 
of a mechanical cantilever. 

We can determine the conditions under which the SC 
is finite. First of all from Eq. (O it follows that the spin- 
orbit coupling must be finite, i.e. k so 7^ 0. Further, from 
Eq. ^ one observes, that when 7 = 0, the Bloch ampli- 



tudes do not depend on j, vSXq j ; i t k B (di t j) 

DVR 



and since [u^ k (d\)] 



,DVR 

n,fc E 



(di), 



(d\) (time-reversal sym- 



metry), the function Fk B ,k' B becomes even with respect to 
its arguments. Then from Eq. ([5]) one gets zero SC. Thus 
the second condition is the presence of the orbit-orbit 
coupling. Finally, since for a symmetric periodic poten- 
tial the Bloch amplitudes are real functions, we conclude 
that the function Fk Bt k B is also real in this case, that is 
ImfJ^jj t f.i ] =0. As a result the third condition is the 
presence of spatial asymmetry. 

Below we present corresponding numerical results. 
All energies and frequencies are given in units of Hloq 
and loq, respectively. The parameters are taken for 
an InGaAs/InP quantum wire: Huq = 0.9 meV; a = 
h 2 k so /m = 9.94- 10 -12 eV-m; m — 0.037mo, respectively. 
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FIG. 3: (Color online) Non-equilibrium spin current, J^ e s> as 
a function of the spin-orbit coupling strength, fc so , for different 
values of the orbit-orbit coupling strength, 7. The driving 
amplitude and viscosity coefficient are F = 2Hcjo/L, rj = 0.5. 
The other parameters are as in Fig. [2] 

For k so L — 7r/2 one gets L = 0.32 /urn. 

The dependence of the NESC on the amplitude 
of the external driving is shown in Fig. [2] for the 
asymmetric periodic potential (see inset) U(x) = 

12n=o V n cos ( 2nnx / L - 4> n ) with V = 4, V\ = -V , 
V 2 = 3.89, 0o = 02 = 0.0, 0i = 1.9. The gap between 
the Bloch bands with I = 1 and I = 2 is AEx2 ~ 10.5. 
In Fig. [U-FL, h£l < AE12 that is the numerical results 
are consistent with the theoretical model assumptions. 
As an example we have used an Ohmic bath with the 
spectral density J(lu) — ?7a;exp(— u>/u> c ), where the vis- 
cosity coefficient (in units of mui ) is 77 = 0.25, 0.5, 0.75, 
and the cutoff frequency is uj c = 10. As it can be seen, 
the NESC has an oscillating nature. However, the oscil- 
lation amplitude goes down when the driving increases. 
Physically such behavior can be attributed to an effective 
renormalization of the band structure in a high-frequency 
electric field [l5j . The group velocity decreases in a non- 
monotonous way which due to RSOI slows down the spin 
kinetics. For increasing values of r\ the dissipation in- 
duced decoherence in the system gets more pronounced. 
The system becomes more classical and thus the tunnel- 
ing processes become less intensive. This leads to the 
spin current reduction which one observes in Fig. [21 

In Fig.[3]the NESC is plotted versus k so L while 7 plays 
the role of a parameter. The oscillations of the NESC 
have minima located at nG/2 where n = 0, 1, 2, . . ., and 
G is the reciprocal lattice vector. Physically this reflects 
the fact that for those values of fc so the Rashba split be- 
comes minimal due to the periodicity of the energy spec- 
trum in the k-space. The magnitude of these oscillations 
decreases with decreasing orbit-orbit coupling, and the 
current vanishes for 7 = 0. 

In summary, we have studied stationary quantum 
transport in a driven dissipative periodic quasi-one- 
dimensional system with Rashba spin-orbit interaction 
and orbit-orbit coupling. The spin ratchet effect has been 
investigated and an analytical expression for the spin cur- 



rent has been derived and analyzed. This analysis has 
revealed that for the case of moderate-to-strong dissipa- 
tion the necessary conditions for non-vanishing spin cur- 
rents are the spatial asymmetry of the periodic potential 
as well as a finite strength of the spin-orbit interaction 
and orbit-orbit coupling. It has been demonstrated that 
in a dissipative system equilibrium spin currents can ex- 
ist. Our numerical calculations have shown characteristic 
oscillations of the spin current as a function of the am- 
plitude of the driving force and the spin-orbit coupling 
strength. Finally, we note, that since the spin current has 
the in-plane polarization, it can be efficiently measured 
by a magneto-optic Kerr microscope using the cleaved 
edge technology as suggested recently in Ref . [2(| . 
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